MEASURING DISLOCATION DENSITY IN ALUMINUM WITH RESONANT 

ULTRASOUND SPECTROSCOPY 
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Dislocations in a material will, when present in enough numbers, change the speed of propaga- 
tion of elastic waves. Consequently, two material samples, differing only in dislocation density, will 
have different elastic constants, a quantity that can be measured using Resonant Ultrasound Spec- 
troscopy. Measurements of this effect on aluminum samples are reported. They compare well with 
the predictions of the theory. 
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I. INTRODUCTION 



II. EFFECTIVE VELOCITY OF ELASTIC 
WAVES IN A DISLOCATION-FILLED MEDIUM 
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A series of papers by Maurel et al. [2004a; 2004b; 
2005a; 2005b; 2006; 2007a; 2007b] have constructed a 
detailed theory of the interaction of elastic waves with 
dislocations in elastic, homogeneous and isotropic, solids. 
This has been done both in two and three dimensions. 
Dislocations have ben considered both in isolation as well 
as in large numbers. In the latter case, generalization of 
the Granato & Lucke [1956a; 1956b; 1966] theory has 
emerged with results for change in wave propagation ve- 
locity and attenuation length that clearly distinguish be- 
tween longitudinal (acoustic) and transverse (shear) po- 
larizations. These results are in satisfactory agreement 
with laboratory measurements of acoustic attenuation 
using Resonant Ultrasound Spectroscopy — RUS [Ledbet- 
ter & Fortunko, 1995; Ogi et al, 1999; 2004; Ohtani et 
al., 2005]. The case of an isolated dislocation in a half- 
space has also been studied, as well as the case of low 
angle grain boundaries mimicked as dislocation arrays. 



A natural development of the above ideas and results 
is to ask whether RUS can be turned into a practical 
tool to measure dislocation densities in materials. In 
order to do this, it is necessary to validate whatever 
results are obtained using RUS with time-tested, but 
more involved, High Resolution Transmission Electron 
Microscopy [Williams & Carter, 2004; Arakawa et al., 
2006; Robertson et al., 2008]. A preliminary step towards 
that aim is to perform RUS measurements on a number 
of samples of a given material, one as received from the 
provider, and the others after cold rolling, or annealing, 
so as to have significantly different dislocation densities. 
This paper reports a simplified form of the theory, as well 
as the first such measurements, using aluminum. 



We study an homogeneous, isotropic, three- 
dimensional, infinite clastic medium of density p, 
whose state is described by a vector field u(x,t), the 
displacement u at time t of a point whose equilibrium 
position is at x. In the absence of dislocations, the 
displacement u obeys the wave equation 



dt 2 



Cijkl 



d 2 u k 
dxjdxi 







(1) 



with Cijki = ^SijSki+p(SikSji+SnSjk) the tensor of elastic 
constants, and i, j, k, ■ ■ ■ = 1,2, 3. A consequence of this 
equation is that the medium allows for the propagation of 
longitudinal (acoustic) and transverse (shear) waves with 
propagation velocity cl = \f{\ + 2^)/p and ct = \fp~Jp-, 
respectively. Their ratio 7 = cl /ct > 1 is always greater 
than one. 

Dislocations are modeled as one dimensional objects 
("strings", [Koehler, 1952; Granato & Lucke, 1956a,b]) 
X(s,t), where s is a Lagrangean parameter that labels 
points along the line, and t is time, of length L, pinned 
at the ends, whose equilibrium position is a straight line. 
They are characterized by a Burgers vector 6, perpen- 
dicular to the equilibrium line. Their unforced motion is 
described by a conventional vibrating string equation 



5 2 A, BXi 9 2 X, 



(2) 



where the mass per unit length m and line tension T are 
given by [Lund, 1988] 
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where S and 5q are external and internal cut-offs. The co- 
efficient B is a phenomenological term that describes the 
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internal losses of the string due to, for example, interac- 
tions with phonons and electrons. We shall only consider 
glide motion, that is, motion parallel to the Burgers vec- 
tor b. 

When elastic waves and dislocations interact, both 
Eqns. ([1]) and ((2]) acquire right hand side — source — 
terms, whose structure has been discussed in detail by 
Maurel et al. [2005a,b] . Elastic waves in the presence of 
N dislocations are best described not in terms of par- 
ticle displacement u but in terms of particle velocity 
v = du/dt and the wave equation (fl]) becomes 



' dt 2 



d 2 v k 
dxjdxi 



Si 



(4) 



the frequency of the fundamental mode of the string with 
fixed ends. 

Maurel et al. [2005b] have provided two derivations 
of effective velocities for elastic waves described by Eqn. 
(j8j. Here we give a third, with a reasoning similar to 
the one used to study waves in plasmas [Stix, 1992]: The 
right-hand-side is smoothed through the replacement 
of the discrete sum over dislocation segments by an in- 
tegral over space with a continuous density n(x) of dis- 
location segments, and the tensor M^ MJ^ by its angular 
average, assuming all directions equally likely. The last 
operation can be found in Appendix C of Maurel et al. 
[2005b]. Eqn. thus becomes, in the case of uniform 
dislocation density n{x) = n, 



where the source term Si is given by 



N f 

Si(x,t) = c ijkl e mnk y2 / ds X^(s,t)T n b, 



x A-S(x-X n (s,t)). 



where e mn k is the completely antisymmetric tensor of or- 
der three, and f is a unit tangent along the dislocation 
line. The string equation @ is written for the compo- 
nent of motion along the glide direction, X = X k t k and, 
loaded by a Peach-Koehler [1950] force it becomes 

mX(s,t)+BX(s,t)-rX"(s,t) = fib M lk d t u k (X,t) : (6) 

with Mi k = tin k + t k ni and n = f x t a unit binormal vec- 
tor. Overdots mean time derivatives, and primes mean 
derivatives with respect to s. 

At this point it becomes profitable to go to the fre- 
quency domain. The loaded string equation ^ can 
be solved in terms of normal modes, and the solution 
plugged into the right hand side of the wave equation 
In the long wavelength limit, A ^ L, and for small 
string displacements, the result of this operation is 
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where 
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8 ^bfS{uj)nL d 2 v k 
2~ T5"( 3 "»fc"*i + q„ a _. ( 10 ) 



it 2 m u> 15 



dxjdx 



(5) In wave number space this is an equation 
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The dispersion relation for the elastic waves in this aver- 
aged medium is given by the vanishing of the determinant 
oiO: 



det O = 



(12) 



For frequencies smaller than the fundamental frequency 
of the string, uj coi, and small damping, (u>B/m) <C 
u) 2 , this leads to the following effective longitudinal (vl) 
and transverse (vt) phase velocities: 
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III. RUS MEASUREMENTS 

Resonant Ultrasound Spectroscopy allows precise mea- 
surements of the elastic constants of a sample indepen- 
dent of its symmetries [Migliori et al., 1993; Leisure & 
Willis 1997]. An homogeneous and isotropic material 
is characterized by two independent elastic constants, A 
and n, or equivalently C\\ — pv L and C44 = pro ji. RUS 
is known to give more precise measurements of C44 , and 
therefore a comparison with the theoretical prediction for 
vt is possible. 
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Parameter 


Sample 1 


Sample 2 


Sample 3 


Sample 4 


Sample 5 


Preparation 


Annealed 400°C/10 hrs 


Annealed 400° C/5 hrs 


Original 


Rolled at 33% 


Rolled at 43% 


di [cm] 


1.701 ±0.001 


1.700 ± 0.001 


1.7004 ± 0.0003 


1.696 ±0.001 


1.7004 ± 0.0005 


d,2 [cm] 


1.0010 ± 0.0004 


0.9992 ± 0.0004 


0.9991 ± 0.0005 


1.000 ±0.001 


1.0005 ± 0.0003 


d 3 [cm] 


5.002 ± 0.001 


5.002 ± 0.001 


5.000 ±0.001 


5.000 ± 0.001 


5.001 ± 0.002 


P [g/cm 3 ] 


2.691 ±0.002 


2.692 ± 0.002 


2.696 ± 0.002 


2.692 ± 0.004 


2.691 ± 0.002 


C*n [10 11 Pa] 


0.868 ± 0.038 


0.880 ± 0.023 


0.884 ± 0.028 


0.888 ± 0.026 


0.883 ± 0.023 


C 44 [10 11 Pa] 


0.2680 ±0.0017 


0.2664 ± 0.0007 


0.2661 ± 0.0008 


0.2647 ±0.0007 


0.2643 ± 0.0007 


v L [m/s] 


5679 ± 124 


5719 ± 76 


5725 ± 89 


5744 ± 84 


5728 ± 75 


v T [m/s] 


3156 ± 10 


3146 ±4 


3142 ± 5 


3136 ±5 


3133 ±4 



TABLE I: Dimensions, density, Cu, C44, vl and vt for the five samples. Columns are ordered for increasing expected density 
dislocation. Absolute errors for Cn and C44 are computed as the standard deviations of a set of 135 = 5 • 3 3 results: for each of 
the 5 positions per sample, the elastic constant seed values used in RUS are varied 3 times each, within an interval close to the 
final measured values (±20%). The density is also varied 3 times within (p) ± 5p, giving a total of 27 possible combinations. 
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FIG. 1: (a) Typical frequency spectrum, showing the first 14 resonances for sample N° 3 of Table I. The first resonant mode 
is at ~ 18 kHz, the last one at ~ 101 kHz. (b) Shear wave velocity for the five samples under study. The arrow indicates the 
direction in which dislocation density is expected to increase. 



We have taken five aluminum samples cut from the 
same bar. One, as bought, one cold-rolled at 33%, 
one cold rolled at 43%, one annealed at 400°C, for 5 
hours, and another for 10 hours. Longer annealing means 
lower dislocation density, and stronger cold-rolling means 
higher dislocation density. The five samples were shaped 
as parallelepipeds, with dimensions as in Table I. Samples 
are labeled 1, 2, 3, 4 and 5 for increasing expected density 
dislocation. A RUS apparatus was built in-house [Caru, 
2007; Jara, 2007], and used to measure the two elastic 
constants of the five samples. The sample-apparatus con- 
tact force is small, of the order of 0.1 N. A typical spec- 
trum is shown in Figure la. Each resonant frequency was 
measured for several ultrasonic driving amplitudes (typi- 
cally five) in order to verify that the resonances are well in 
the linear acoustic regime. Additionally, each sample was 
placed five times in the apparatus in order to reduce er- 
rors due to slight dependence of the resonant frequencies 
on the contact load and positioning with respect to the 
ultrasonic receiver. The measured elastic constants Cn 



and C44 as well as the shear and longitudinal wave veloc- 
ities are also given in Table I. No clear tendency in vl is 
observed. Within experimental errors, it is almost con- 
stant. However, vt shows a clear decreasing tendency. 
It is plotted in Figure lb versus sample number. The 
difference between the errors of both wave velocities is 
consistent with the fact that RUS is much more precise 
for C44, basically because resonant frequencies depend 
strongly on C44 and weakly on Cn. 

IV. DISCUSSION AND CONCLUSIONS 

For simplicity we shall assume, a common assumption, 
A = /i, so that 7 2 = 3. Independent measurements for 
the original sample give 7 2 = 3.3, and this ratio does not 
change significantly from sample to sample. 

Using © we have 

/x& 2 47T 7 2 _ x / 5 \ 

- = y-T ln UJ 
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so there is only a dependence on the ratio of cut-offs. 
Taking S = 10 6 5o we get the following expression for the 
fractional change in shear wave velocity: 

vt-ct _ 16 3 nL 3 

c T 57r 3 2 61nl0' 1 ' 

The data of Figure lb are consistent with 

nL 3 ~ 0.3 (16) 

and with a (linear) trend in agreement with theory 
namely that the higher the dislocation density the lower 
the effective speed of shear waves. 

Taking L ~ 100 nm as a typical dislocation length 
would give a variation in dislocation density among the 
various samples of order nL ~ 10 9 mm~ 2 , a conclusion 
that it should be possible to test by direct measurement 
with High Resolution Transmission Electron Microscopy 
(HRTEM). Work along this direction is in progress. 
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